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Abstract. We study the spectral properties of a low angular momentum flow 
as a function of the shock strength, compression ratio, accretion rate and flow 
geometry. In the absence of a satisfactory description of magnetic fields inside 
the advective disk, we consider the presence of only stochastic fields and use the 
ratio of the field energy to the gravitational energy density as a parameter. We not 
only include 'conventional' synchrotron emission and Comptonization by Maxwell- 
Bolzmann electrons in the gas, but we also compute these effects due to power-law 
electrons. For strong shocks, a bump is produced due to the post-shock flow. A 
power-law spectral components due to the thermal and non-thermal electrons appear 
after this bump. 
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1. Introduction 

In Chakrabarti & Titarchuk (1995, hereafter CT95), a two component 
advective flow (TCAF) solution was introduced (see, Chakrabarti, this 
volume). The sub-KepIerian flow has a steady shock but depending on 
the flow parameters it could be oscillating or propagating. The post- 
shock region is called CENtrifugal pressure supported BOundary Layer 
or CENBOL. Given that the accretion flow largely passes through this 
shock, it is likely that the hot electrons may be accelerated by the shock 
just as high energy cosmic rays are produced by the transient super- 
novae shocks (Bell 1978ab, Longair 1981). We assume this process in 
computing the spectra and show that the presence of shocks make 
their mark distinctly. Our solution also produces gamma-rays naturally 
since shocks are generally produced due to the centrifugal barrier (e.g., 
Chakrabarti, 1996). It is possible that occasionally an additional source 
of high energy photons, perhaps in the form of shocks in outflows and 
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jets, is necessary to explain the complete spectra. Details of our model 
are presented in Mandal &; Chakrabarti (2004). 



2. Relevant heating and cooling processes and the solution 

procedure 

Protons loss energy through Coulomb interaction and inverse bremss 
-trahlung. Electrons are heated through this Coulomb coupling. Elec- 
trons are cooled by bremsstrahlung, cyclo-synchrotron and Comptoniza- 
tion of the soft photons due to cyclo-synchrotron radiation. The Comp- 
tonization is computed by using this cooling term augmented by the 
enhancement factor. In the present situation, the flow passes through 
the shocks and thus the electrons will have a power-law distribution 
Thus, both the thermal and the non-thermal electrons are present. 

The geometry of the flow is chosen to be conical for simplicity and 
the flow surface makes an angle with the z-axis. We assumed the form 
of infall velocity v{x) = vqx~^^'^. At the shock, x = Xs, the velocity is 
reduced by a factor R, the compression ratio, i.e., V-^-{xs) = V-{xs)/R. 
The number density also goes up by a factor R. 

As for the initial conditions, we flx the outer boundary at a large 
distance (say, lO^rg) and supply matter (both electrons and protons) 
with the same temperature (say, Tp = Te = lO^K). By integrating en- 
ergy equation, wc obtained the temperature. The radiation emitted by 
the flow through bremsstrahlung and synchrotron radiation are inverse 
Comptonized by the hot electrons in the flow. 

A shock of compression ratio R causes the formation of power-law 
electrons of slope (Bell, 1978ab), p = [R + '^)/{R — 1). This power- 
law electrons produce a power-law synchrotron emission with an index 
(Longair, 1981) q = {1 —p)/2. The power-law electrons have energy 
minimum at i^niin and energy maximum at £^max- Smux is obtained self- 
consistently by conserving the number of power-law electrons and by 
computing the number of scatterings that the electrons undergo inside 
the disk before they escape and f^min is obtained from the temperature 
of the injected electrons. This temperature is obtained self-consistently 
through our integration procedure. 

In presence of magnetic field (assumed to be the ratio of the mag- 
netic energy density to gravitational energy density), the pre-shock 
flow and the CENBOL would emit synchrotron radiation. We assume 
that the number of black body photons generated by the electrons 
obeying power-law distribution (assumed to be of constant fraction 
C of thermal electrons) below self-absorption frequency is the same 
fraction ( of the total black body photons. We followed the procedures 
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Figure 1, (left) Wedge shaped flow model assumed in this paper which included 
a shock at the inner region, (right) Variation of proton and electron tempera- 
tures (marked as Tp and Te) as functions of dimensionless radial distance x for 
m = 0.00059. Other parameters are M = WMq, R = 3.9, 6 = 77°, C = 0.7 and 
Kg = 10. Closer to the black hole, the electrons cool more rapidly than the protons. 

presented in CT95 and Titarchuk &; Lyubarskij (1995) while computing 
the Comptonized spectrum due to scattering from thermal electrons 
and power-law electrons. However, unlike CT95 where a Keplerian 
disk was supplying soft photons, here the source is everywhere, i.e., 
distributed throughout. This has been taken into account. At the end, 
we add contributions from all the components to get the net photon 
emission from the flow. 

In the next Section, we present the spectral characteristics by vary- 
ing parameters such as 0, Xg, R, C and m. We use a black hole of mass 
1OM0. 

3. Results and interpretations 

In Fig. 1 we show our simplified flow geometry to capture the essential 
physics and show how the temperature of the electrons and protons 
differ as the flow moves closer to the black hole. The deviation is 
more prominent in the post-shock region. Here, the radial distance x is 
measured in units of the Schwarzschild radius Vg. 

In Fig. 2, we show a typical spectrum obtained from our solution. 
The curve marked 1 is due to pre-shock flow and the curve marked 2 is 
due to the post-shock flow. The synchrotron self-absorption frequency 
at fa is marked. The bump due to post-shock flow is at frequency 

t'bump ~ t'injil + I^TTS-]""'" , where, finj is the frequency of the 
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Figure 2. A typical complete spectrum from a sub-Keplerian accretion flow with all 
the contributions shown for the same parameters as in Fig. lb. The curve marked '1' 
is due to Comptonization of the synchrotron emission from the pre-shock region. The 
curve marked '2' is due to the Comptonized spectra from the post-shock (CENBOL) 
region. The curve marked '3' indicates the net radiation from the flow. Vb. represents 
synchrotron self-absorption frequency for the CENBOL. 
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Figure 3. Variation of the total emitted spectrum when the shock location Xs is 
varied. The solid, dotted and dashed curves are drawn for Xe = 10, 40, 80 respec- 
tively. The spectrum becomes softer when Xa is reduced. Here, rh = 0.1 has been 
chosen. 



dominant photons from the pre-shock flow, R is the compression ratio 
of the shock located at Xg. The power law is produced by the power- law 
electrons from the shock acceleration. 

In Fig. 3, we compare three cases with shock locations at Xg = 10, 40 
and 80 respectively. The shock locations can vary with specific angu- 
lar momentum in the flow which in turn is governed by the viscosity 
parameter (Chakrabarti, 1990; 1996). 

With the rise in the accretion rate of the sub-Keplerian disk, the 
spectrum is likely to become harder. This can be seen more clearly in 
Fig. 4, where we vary the accretion rate. A higher accretion rate is 
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Figure 4- Variation of the spectra with the dimensionless accretion rate. Solid, dot- 
ted and dashed curved are for rh = 0.05, 0.5 and 1.5 respectively. The spectrum 
becomes harder when the accretion rate is increased. 




Figure 5. Example of a typical soft state and a hard state spectra of a galactic 
black hole candidate os mass IOM0. The parameters for the hard state are: m — 2.4, 
R = 4, Xs = 250, Q = 80, C = 0.24. The parameters for the soft state are: rh = 0.02, 
R = 4,x^ = 100, e = 80, C = 0.3. 



expected to increase the density of the flow and it becomes difficult to 
cool the matter by the soft photons of the synchrotron radiation. The 
Comptonized spectrum becomes harder. At lower accretion rate, it is 
easy to cool the flow and the spectrum becomes softer. The curves are 
drawn for: rh = 0.05 (solid), 0.5 (dotted) and 1.5 (dashed). 



shockidenti.tex; 2/02/2008; 15:48; p. 5 



6 



S. Mandal and S.K. Chakrabarti 



4. A typical soft/hard spectra in gamma-ray region 

Typical broad band soft state and hard state spectra are presented in 
Fig. 4 of Chakrabarti (this volume). In Fig. 5 of Chakrabarti (this vol- 
ume) behaviour at high energy has been shown for Cyg X-1, GROJ1719- 
24 and GROJ0422+32 (taken from Ling k Whcaton, 2004). In Fig. 5 
we present the photon flux in soft and hard states as obtained by our 
solution with shocks. No jet contribution was added. Not surprisingly, 
they look similar. However, a second kind of high 7 intensity state 
was reported (sec, Case et al. 2004) in which non-thermal power-law 
has been observed from a few tens of keV to few MeVs (see. Fig. 4 
of Chakrabarti, this volume). It has been difficult to obtain such a 
behaviour without invoking jets or other mechanisms. 
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